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Comparisons were made between the life prediction formulas of Lundberg and Palmgren, Ioannides and Harris, and Zaretsky and full- 
scale ball and roller bearing life data. The effect of Weibull slope on bearing life prediction was determined. Life factors are proposed 
to adjust the respective life formulas to the normalized statistical life distribution of each bearing type. The Lundberg- Palmgren method 
resulted in the most conservative life predictions compared to Ioannides and Harris and Zaretsky methods which produced statistically 
similar results. Roller profile can have significant effects on bearing life prediction results. Roller edge loading can reduce life by as 
much as 98 percent. The resultant predicted life not only depends on the life equation used but on the Weibull slope assumed, the least 
variation occurring with the Zaretsky equation. The load-life exponent p of 10/3 used in the ANSI/ABMA/ISO standards is inconsis- 
tent with the majority roller bearings designed and used today. 

Keywords: Rolling-Element Bearings; Life Prediction Methods; Rolling-Element Fatigue 


INTRODUCTION 

Rolling bearing technology has evolved over 4000 years to 
the present. H.T. Morton [ 11 in his 1965 book, “Anti-Friction Bear- 
ings/’ describes the evolution of rolling bearing technology. 
Although rolling bearing technology continued to develop 
throughout the first half of the 19th century, it is not until the 
invention of the pedal bicycle that the rolling-element bearing 
industry becomes established. In 1868, A.C. Cowper made a 
bicycle with ball bearings. According to Morton, W. Bown of 
Coventry, England, appears to be the most successful bearing 
manufacturer. In 1 880, he had a contract to produce 12 ball bear- 
ings a day for Singer and Company, a bicycle manufacturer. 
By 1920, most types of rolling-element bearings used today were 
in production. 

By the close of the 19th century, the bearing industry began 
to focus on the sizing of bearings for specific applications and 


NOMENCLATURE 

A material-life factor 

C dynamic load capacity, N (lbf) 

c critical shear stress-life exponent 

e Weibull slope 

F probability of failure, fraction or percent 

j{X) probability of survival function 

h exponent 

L life, number of stress cycles or hr 

L x 0 10-percent life or life at which 90 percent of a 

population survives, number of stress cycles or hr 
N life, number of stress cycles 

n maximum Hertz stress-life exponent or number of 

components, elemental volumes 
P normal or equivalent radial load, N, (lbf) 

p load-life exponent 


determining bearing life and reliability. In 1896, R. Stribeck [2] 
began fatigue testing full-scale bearings. In 1912, J. Goodman 
[3] published formulas based on fatigue data that would compute 
safe loads on ball and cylindrical roller bearings. 

The most influential person in developing life prediction meth- 
ods for ball and roller bearings was A. Palmgren [4] in Sweden. 
His work with G. Lundberg [5,6] published in 1947 and 1952, 
modified and developed the theoretical basis for Palmgren’s 1945 
formulas [7] by using the Weibull distribution function [8,9] . Their 
work resulted in the International organization for Standardiza- 
tion (ISO) and the American National Standards Institute (ANSI)/ 
Anti-Friction Bearing Manufacturers Association (AFBMA) (now 
American Bearing Manufacturers Association (ABMA)) stan- 
dards for the load ratings and life of rolling-element bearings 
[10 to 12]. The basis for the Lundberg -Palmgren theory and the 
ISO-ANSI/ABMA standards were outlined in a German paper 


S probability of survival, fraction or percent 

S max maximum Hertz stress, GPa (ksi) 

V stressed volume, m 3 (in. 3 ) 

X load, time, or stress 

Z depth to maximum critical shear stress, m (in.) 

a stress or strength, GPa (ksi) 

r critical shear stress, GPa (ksi) 

T u fatigue limit, GPa (ksi) 

Subscripts 

act. actual life based on laboratory or field data 

n number of components or elemental volumes 

pred. predicted life based on designated life prediction method 
P designates characteristic life or stress 
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presented by Palmgren [13] in 1924. This paper is the first to 
propose the L ]0 life as the basis for life prediction and also the 
“linear damage rule” to determine the effect of varying load. 

As rolling bearing design and steels improved, the ISO stan- 
dard and the Lundberg- Palmgren equations were found to under 
predict bearing lives. In order to account for these differences in 
life, ASME published “Life Adjustment Factors” in 1971 [14], 
These life adjustment factors were modified and expanded by 
the STLE in 1991 [15]. In 1985, E. Ioannides andT. Harris [16] 
proposed a modification of the Lundberg- Palmgren theory based 
upon a material fatigue limit and discrete finite elements. The 
concept of a fatigue limit was first introduced by Palmgren in 
1924 [13] and then subsequently disregarded by him. 

In 1987, E.V. Zaretsky [17] proposed a Weibull-based life 
theory using discrete finite elements that accounts for deviations 
in the accepted material Hertz stress-life relation. Zaretsky, 
together with J. V. Poplawski and S. Peters, applied his life theory 
to ball and roller bearings [ 1 8, 19]. 

Other researchers have proposed modifications of the 
Lundberg-Palmgren theory based upon both a fatigue limit and 
debris contamination [20 to 23]. T.E. Tallian [24 to 26] has 
developed a mathematical model also based on Lundberg- 
Palmgren but which is fitted to a bearing data base. 

The work reported and discussed herein is based in part on our 
work presented in [18,19]. We summarize and compare the life 
prediction formulas of Lundberg and Palmgren, Ioannides and 
Harris, and Zaretsky to each other and the ISO standards as well 
as with full-scale ball- and roller-bearing life data. 

LIFE THEORIES 
Weibull Equation 

Fracture Strength — In 1939 W. Weibull [8,9] published 
two papers that describe a statistical approach to determine the 
strength of solids. Weibull postulated that the dispersion in mate- 
rial strength for a homogeneous group of test specimens could be 
expressed according to the following relation: 

lnln(l/5) = ^ln[x/X p ] (1) 

where X = o and Xp = Gp . 

Equation (2) relates specimen survival S to the fracture (or rup- 
ture) strength <r. When In ln( 1/5) is used as the ordinate and In cr 
as the abscissa and fracture (and fatigue) data are assumed to plot 
as a straight line. The slope (tangent) of this line is referred to as 
the Weibull slope or Weibull modulus usually designated by the 
letter e or m. The plot itself is referred to as a Weibull plot. 

By using a Weibull plot, it becomes possible to estimate a 
cumulative distribution of an infinite population from an extremely 
small sample size. The Weibull slope is indicative of the dispersion 
of the data and its density (statistical) distribution. Weibull slopes 
of 1, 2, and 3.57 are indicative of exponential, Rayleigh, and nor- 
mal (Gaussian) distributions, respectively [8]. 


The scatter in the data is inversely proportional to the Weibull 
slope, that is, the lower the value of the Weibull slope, the larger 
the scatter in the data and vice versa. The Weibull slope is also 
liable to statistical variation depending on the sample size (data 
base) making up the distribution [27]. The smaller the sample 
size the greater the statistical variation in the slope. 

Weibull [8,9] related the material strength to the volume of the 
material subjected to stress. If we imagine the solid to be divided 
in an arbitrary manner into n volume elements, the probability of 
survival for the entire solid can be obtained by multiplying the 
individual survivabilities together as follows: 

5=5 r S 2 -5 3 ...5 w (2) 

where the probability of failure F is 

F = l-S (3) 

Weibull [8,9] further related the probability of survival S, the 
material strength a, and the stressed volume V, according to the 
following relation: 

lni = J y /(XMV (4) 

where 

f(X) = c e (5) 

For a given probability of survival 5, 



From Eq. (6) for the same probability of survival the compo- 
nents with the larger stressed volume will have the lower strength 
(or shorter life). 

Fatigue Life — In conversations with E.V. Zaretsky on 
January 22, 1964, W. Weibull related how he had suggested to 
his contemporaries A. Palmgren and G. Lundberg in Gothenburg, 
Sweden, to use his equation to predict bearing (fatigue) life where 

/(X) = x c N e (7) 

and where r is the critical shear stress and N is the number of 
stress cycles to failure. 

In the past we have credited this relation to Weibull. How- 
ever, there is no documentation of the above nor any publication 
to the authors’ knowledge of the application of Eq. (7) by Weibull 
in the open literature. However, in [18] we did apply Eq. (7) to 
Eq. (4) where 
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The parameter c/e is the stress -life exponent. This implies that and for point contact 

the inverse relation of life with stress is a function of the life 
scatter or data dispersion. 

From Hertz theory, V and x can be expressed as a function 
of S max 1 1 8] and substituting L for N 


(14b) 


L = A\ 


r i il/e 


V 

c/e 

rii 

^ T_ 


LvJ 


S n 

J max 


From Lundberg and Palmgren [5], using the values of c and e 
previously discussed and h = 2.33, then from Eqs. ( 14a) and ( 14b), 
n equals 8.1 and 9 for line and point contact, respectively. 


From [18], solving for the value of the exponent n for line con- 
tact (roller on raceway) from Eq. (9) gives 

c + 1 

n = (10a) 

e 


and for point contact 


Ioannides-Harris Equation 

Ioannides and Harris [16], using Weibull [8,9] and Lundberg 
and Palmgren [5,6] introduced a fatigue-limiting stress where from 
Eq. (4) 


f(X) = 


(x-xj c ^ g 

Z h 


(15) 


c + 2 

n = (10b) 

e 

Using the value of e = 1.11 and c/e = 9.3 from Lundberg and 
Palmgren [5], n equals 10.2 and 11.1 or line and point contact, 
respectively. 


Lundberg-Palmgren Equation 

In 1947, G. Lundberg and A. Palmgren [5] applied Weibull 
analysis to the prediction of rolling -element bearing fatigue life. 
The Lundberg-Palmgren theory expressed fiX) in Eq. (4) as 

T c N e 

/w= ^r (id 

where ris the critical shear stress, N is the number of stress cycles 
to failure, and Z is the depth to the maximum critical shear stress 
in a concentrated (Hertzian) contact From Eqs. (4) and (11) 


From Eqs. (4) and (15) 


N 


\c/e 




1/e 


[Z] 


hie 


(16) 


Equation (16) is identical to that of Lundberg and Palmgren 
(Eq. (12)) except for the introduction of a fatigue-limiting stress. 
Equation (16) can be expressed as a function of 5 m where 



Ioannides and Harris [ 1 6] use the same values of Lundberg and 
Palmgren for e , c, and h . If x u equal 0, then the values of the 
exponent n are identical to those of Lundberg and Palmgren (Eqs. 
(14a) and (14b)). However, for values of x u > 0, n is also a func- 
tion of(r- x u ). 


N ~ 


“r 

de 

rii 

_T_ 


LvJ 


{Z\ Me 


(12) 


From Hertz theory, V, r, and Z can be expressed as a function of 
^max’ substituting L for N 


/ * \c/e / | xl/e* 

L = A - - [Z] 

U ) \v) 1 J 


\hfe 1 


C« 

‘-’max 


(13) 


Solving for the exponent n for line contact gives 

c + l-/z 


n ~ ■ 


(14a) 


Zaretsky Equation 

Both the Weibull and Lundberg-Palmgren equations above 
relate the critical shear stress-life exponent c to the Weibull slope 
e. The parameter c/e thus becomes, in essence, the effective criti- 
cal shear stress-life exponent, implying that the critical shear 
stress-life exponent depends on bearing life scatter or dispersion 
of the data. A search of the literature for a wide variety of mate- 
rials and for nonrolling-element fatigue reveals that most stress- 
life exponents vary from 6 to 12. The exponent appears to be 
independent of scatter or dispersion in the data. Hence, Zaretsky 
[28] has rewritten the Weibull equation to reflect that observation 
by making the exponent c independent of the Weibull slope e, 
where 


f(X) = x ce N e (18) 
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From Eqs. (4) and (18) 



(19) 


Equation (19) differs from the Weibull Eq. (8) and the 
Lundberg-Palmgren Eq. ( 1 2) in the exponent of the critical stress 
r. Zaretsky assumes based upon experience that the value of the 
stress-exponent c- 9. 

Lundberg and Palmgren (5) assumed that once initiated, the 
time a crack takes to propagate to the surface and form a 
fatigue spall is a function of the depth Z to the critical shear stress. 
Hence, by implication, bearing fatigue life is crack propagation 
time dependent. However, rolling-element fatigue life can be 
categorized as “high-cycle fatigue.” (It should be noted that at 
the time (1947) Lundberg and Palmgren published their theory, 
the concepts of “high cycle” and “low cycle” fatigue were only 
then beginning to be formulated.) Crack propagation is an 
extremely small time fraction for the total life or running time of 
the bearing. The Lundberg-Palmgren relation implies that the 
opposite is true. To decouple the dependence of bearing life and 
crack propagation rate, Zaretsky [8] dispensed with the Lundberg- 
Palmgren relation of L ~ Z hie in Eq. (19). 

Equation (19) can be written as 



( 20 ) 


From Ref. [ 18], solving for the value of the exponent n, for line 
contact from Eq. (20) gives 


1 

n = c + — 
e 

and for point contact 


(21a) 



Life factor, LAjso 
( b) Cylindrical roller bearings. 

Figure 1 . — Comparison of life theories for rolling-element 
bearings with ISO standard using closed form solution. 
Fatigue limiting shear stress for 145 assumed for 
loannides-Harris, 276 MPa (40 ksi). 

similar to that of Zaretsky if the exponents are chosen to be iden- 
tical. Comparison of the Lundberg-Palmgren equation and the 
Zaretsky equation with the ANSI/ABMA/ISO standards is shown 
in Fig. 1 for ball and cylindrical roller bearings. The theoretical 
lives were normalized to a maximum Hertz stress of 4. 14 GPa 
(600 ksi) and subsequently normalized to the calculated ANSI/ 
ABMA/ISO standards at each stress level. For the Ionanides- 
Harris comparison shown in Fig. 1, a fatigue limiting stress of 
276 MPa (40 ksi) was assumed. For ball bearings, the ANSI/ 
ABMA/ISO standard and the Lundberg-Palmgren equation give 
identical results. For roller bearings, the results are not identical. 

From Lundberg and Palmgren [5], the L [0 life of a bearing 
can be determined from the equation: 

ko =[C/P] p (22) 


n = c + — (21b) 

e 

where c = 9 and e = LI 1, n - 9.9 for line contact and n = 10.8 
for point contact. 

Although Zaretsky [28] does not propose a fatigue-limiting 
stress, he does not exclude that concept either. However, his 
approach is entirely different from that of Ioannides and Harris 
[16]. For critical stresses less than the fatigue-limiting stress, the 
life for the elemental stressed volume is assumed to be infinite. 
Thus, the stressed volume of the component would be affected 
where L = VV^ e . As an example, a reduction in stressed volume 
of 50 percent results in an increase in life by a factor of 1 .9. 

RESULTS AND DISCUSSION 
Life Theory Comparison 

The loannides-Harris analysis without a fatigue limit is iden- 
tical to Lundberg-Palmgren analysis. The Weibull analysis is 


where L 10 is the bearing life in millions of race revolutions, C 
is the dynamic load capacity of the bearing or the theoretical load 
that will produce a life of 1 -million race revolutions with a 
90 percent probability of survival, P is the applied equivalent 
radial load and p is the load-life exponent. For point contact, 
p equals n!3. For line contact, p equals rc/2. 

The ANSI/ABMA/ISO standards use a load-life exponent 
p of 10/3 (3.33) for line contact. This results in a value of n equal 
to 6.6 and can account in part for the lower life predictions than 
that experienced in the field. 

From Lundberg and Palmgren [5], the load-life exponent 
p for line contact should be 4. This results in a maximum Hertz 
stress-life exponent n of 8.1 for line contact (see Eq. ( 14a) i. This 
value of n while low is consistent with available but limited data 
[18,29]. 

Lundberg and Palmgren’s justification for a p of 10/3 was 
that a roller bearing can experience “mixed contact,” that is, one 
raceway can experience line contact and the other raceway point 
contact [6]. This may be true in limited roller bearing designs but 
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TABLE 1.— MAXIMUM HERTZ STRESS-LIFE EXPONENT AS 
FUNCTION OF WEIBULL SLOPE FOR FOUR LIFE 
EQUATIONS 1 (From [19].) 


Equation 

Weibull 

slope 

Stress-life exponent, 
n 

Load-life exponent, 

P ... . 

Line 

contact 

Point 

contact 

Line 

contact 

Point 

contact 

ANSI/ 

ABMA/ISO 

1.11 

6.6 

9 

3.33 

3 

Weibull 
eq. (9) 

1.11 

10.2 

11.1 

5.1 

3.7 

1.5 

7.5 

8.2 

3.8 

2.7 

2.0 

5.7 

6.2 

2.9 

2.1 

Lundberg- 
Palmgren, 
eq. (13) 

1.11 

8.1 

9 

4.1 

3 

1.5 

6.0 

6.7 

3 

2.2 

2.0 

4.5 

5.0 

2.3 

1.7 

Ioannides- 
Harris, 
eq. (17) 

1.11 

8.1 

9 

4.1 

3 

1.5 

6.0 

6.7 

3 

2.2 

2.0 

4.5 

5.0 

2.3 

1.7 

Zaretsky, 
eq. (20) 

111 

9.9 

10.8 

5 

3.6 

1.5 

9.7 

10.3 

4.9 

3.4 

2.0 

9.5 

10.0 

4.8 

3.3 


‘No fatigue limit assumed, u equal 0. 


it is certainly not consistent with the vast majority of cylindrical 
roller and tapered roller bearings designed and used today. 

For line and point contact, the maximum Hertz stress-life and 
load-life exponents were determined for the Weibull, Lundberg- 
Palmgren, Ioannides-Harris, and Zaretsky equations as a 
function of the Weibull slope or Weibull modulus. These results 
are summarized in Table L Both the load-life and stress-life rela- 
tions are based upon the value of the Weibull slope which for 
rolling-element bearings is assumed to be 1.11. For Lundberg 
and Palmgren this assumption resulted in their analysis matching 
preexisting life equations of Palmgren [7] and their nonpublished 
bearing life data base. 

As shown in Table 1, both the load-life and stress-life 
relations of Weibull, Lundberg and Palmgren, and Ioannides and 
Harris reflect a strong dependence on the Weibull slope. The 
existing rolling -element fatigue data reported by Parker and 
Zaretsky [18] reflect slopes in the range of 1 to 2 and some cases 
higher or lower. If the slope were factored into the equations then, 
as shown in Table 1, the stress-life (load-life) exponent signifi- 
cantly decreases with increases in Weibull slope whereby the 
relation no longer matches reality. 

The Zaretsky equation that decouples the dependence of the 
critical shear stress-life relation and the Weibull slope shows only 
a slight variation of the maximum Hertz stress-life exponent n 
and Weibull slope. The value of n varies between 9.5 and 9.9 for 
line contact and 10 and 10.8 for point contact for Weibull slopes 
between 2 and 1.11. 


These results would indicate that for a ninth power Hertz 
stress-life exponent for ball bearings, the Lundberg-Palmgren 
formula best predicts life. However, for a 12th power relation 
reflected by modem bearing steels, the Zaretsky formula based 
on the Weibull equation is superior. 

Poplawski, Zaretsky, and Peters [19] analyzed the effect of 
roller profile on cylindrical bearing life. The results of their 
analysis are summarized in Table 2 for the Lundberg-Palmgren 
method and the Zaretsky method. 

With a closed form solution and not considering edge or stress 
concentrations, the flat roller profile has the longest predicted 
life followed by the end-tapered profile, the aerospace profile 
and the crowned profile, respectively. The full-crowned profile 
produces the lowest lives. While there are life differences 
between the end- tapered profile and the aerospace profile, these 
differences may not be significant. The effect of edge loading on 
the flat roller profile is to reduce life at higher loads by as much 
as 98 and 82 percent at lower loads. The actual percentage calcu- 
lated depends on the analysis used. The effect of roller profile 
needs to be incorporated into the life prediction methodology. 

Comparison With Bearing Life Data 

Harris [30,3 1] analyzed 62 rolling -element bearing endurance 
sets. These data were obtained from four bearing manufacturers, 
two helicopter manufacturers, three aircraft engine manufactur- 
ers, and U.S. Government agency -sponsored technical reports. 
The data sets comprised deep-groove radial ball bearings, 
angular-contact ball bearings, tapered roller bearings and spheri- 
cal roller bearings totaling 7935 bearings. Data for the tapered 
roller bearings were not included by Harris [30,31] because of 
insufficient bearing internal geometry for analysis. Also, data for 
spherical roller bearings and some cylindrical roller bearings were 
not included because the operating conditions were too complex 
for accurate analysis. Of the 62 data sets 51 were analyzed. 
We applied to these 5 1 data sets the Ioannides-Harris life method 
(Eq. (23)) and the Zaretsky life method (Eq. (26)). We used the 
results from Harris [30] for the Lundberg-Palmgren method 
(Eq. (19)) and the modified Harris method without any changes. 
Using Weibull statistics we evaluated these data for each bearing 
type at the L 10 life level and compared the results in Table 3. 

Table 3 compares both the mean and median values of the 
ratio of the actual L 10 life to that predicted for angular-contact 
ball bearings, deep-groove ball bearings and cylindrical roller 
bearings. Ideally, a life ratio (L 10(acl / Lj o<pred >) °f 1 is desirable. 
This indicates that the actual L 10 life (act.) is that which was 


TABLE 2.— EFFECT OF ROLLER PROFILE ON RELATIVE LIFE 1 (From [19].) 
[Weibull slope, 1.11, critical shear stress, 45 .] 


Maximum 
Hertz stress, 
GPa, (ksi) 

ANSI/ 

ABMA/ISO 

standard 

Closed form solution without edge loading i 

Lundberg-Palmgren (eq. (13)) 

Zaretsky (eq. (20)) 

Hat 

End tapered 

Aerospace 


Hat _ 



Crowned 

K sm 

MEM 

■tta 

mm 

■wra 

mm 

isa 


BBtim 

~n 

ma 


ma 


wra 

El 

■wa 

mm 


MEM 

wxm 


Em 

ESI 

eh 

mm 

Em 

ESI 

8.1 

■Egg 

MSfeM 

Ea 


KM 

Em 

eh 

Em 



EE9 

mm 

Em 

in 

sm 

EH 

Em 

Em 

eh 

mm 

WSM 

mm 


El 

Em 

sm 

Em 

2.4 (350) 

• 

- 

2.2 

- 

1.6 


in 


Em 

- 



■nca 


im 


Em 



1 IV iUUAUiiuiu 11CU4. “T. VJI a (UW KM, 

formalized to maximum Hertz stress of 2.4 GPa (350 ksi). 
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TABLE 3.— STATISTICAL COMPARISON OF FIELD AND 
LABORATORY BEARING LIFE DATA TO THAT PREDICTED 


Bearing type 


Weibull 

modulus 

Probability 
of over 
predicting 
life, 
percent 

■ 

Median 

Lundberg-Palmgren method 4 I 

Angular-contact ball brgs 1 

■E9 

6.4 


18 

■ vm i3m >Li i 

Bill 

3 

1.41 

14 

rwrn7iwar?iBT?n— 

20.1 

14.1 

0.98 

5 

Ioannides-Harris method 5 I 

Angular-contact ball brgs 1 

5.7 

2.7 

0.70 

29 

Deep-groove ball brgs 2 


3.3 

1.14 

16 

Cylindrical roller brgs 3 

12.4 

7.8 

0.89 

10 

Modified Ioannides-Harris method 6 | 


4.6 

1.9 

0.65 

36 


■BM 

2.4 

1.55 

17 


mm 

3.2 

0.98 

20 

Zaretsky method 7 

Angular-contact ball brgs 1 

3.1 

1.6 

0.75 

39 

Deep-groove ball brgs 2 

Kfil 

2.2 

1.36 

21 



5.0 

1.02 
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1 14 bearing data sets. 
z 26 bearing data sets 
3 1 1 bearing data sets 

4 From Harris [30]. Contains STLE life factors [15]. 

5 From Eq. (17). Contains STLE life factors [15]. Does not contain 
life factors for material and processing. 

6 From Harris [30]. Contains life factors for lubricant contamination 
only. 

7 From Eq. (20). Contains STLE life factors [15], Normalized to a 
maximum Hertz stress of 4.14 GPa (600 ksi). 

predicted (pred.). Ratios greater than 1 mean that the predicted 
values were lower than the actual values. For life ratios less than 1, 
the predicted values are higher than the actual values. A 
median ratio of 1 would mean that one-half of the predictions 
were high and the other one-half were low. The least incidence of 
over prediction occurs with the Lundberg-Palmgren method. For 
both the Ioannides-Harris method and the Zaretsky method, the 
incidence of over prediction is statistically equivalent. 

The Weibull modulus indicates the dispersion in the 

L 10(act / L 10(pred.) data - ^ hi S her *** value ° f ^ WeibuU modu " 
lus, the less scatter in the life ratios. The lower the value, the 

higher the scatter. For each bearing type, the dispersion in the 
life predictions approximates an exponential distribution. The 
most desirable result would be one where the Weibull modulus 
approaches infinity and the median life ratio is 1 . 

It would be desirable to have the Lio<actY^iO(pred.) rados 
of all bearing types greater than 1 at reasonable levels of 
engineering certainty. In order to accomplish this, additional life 
adjustment factors need to be configured into the bearing life 
calculations. Based upon a statistical analysis using the Weibull 
parameters of Table 3, life adjustment factors were determined 
for each life method and bearing type whereby the predicted 
L l0 life would be equal to or less than the actual L 10 life at a 
designated confidence level. These statically determined life 
adjustment factors are shown in Table 4. 

Table 4 is based upon the bearing life data of Harris [30,31]. 
While it is probable that all design and operating parameters 
necessary to more accurately calculate bearing life were not 


available to Harris [30,31], these bearing life data are the best 
compilation in the open literature. The use of these life factors 
will result in very conservative values of predicted life (L 10< ^ 

A calculated L 10 life multiplied by a life adjustment factor at a 
confidence level of 90 percent should result in an actual life 
( L l(Xact.)) win ** S reater thm 1116 predicted life (L l(Kpn±) ) 
90 out of 100 times regardless of whose life prediction formula 
is used. At a confidence level of 99 percent, when using the life 
adjustment factors, the ratio of ^"KX act. /^"10( pred.) S 1 should 
occur 99 out of 100 times. Whether the designers and users of 
bearings should use these life adjustment factors or variants of 
them needs careful consideration and deliberation. That is for the 
future to decide. 

SUMMARY OF RESULTS 

Comparisons were made between the life prediction formu- 
las for Lundberg and Palmgren, Ioannides and Harris, and 
Zaretsky and their application to full-scale ball and roller bearing 
life data. The effect of Weibull slope on bearing life prediction 
was determined. Life factors are proposed to adjust the respec- 
tive life formulas in order to normalize them to the data base of 
each bearing type. The following results were obtained: 

1 . For the three bearing types studied, the Lundberg-Palmgren 
method resulted in the most conservative life predictions 
compared to that of Ioannides-Harris and Zaretsky. The methods 
of Ioannides-Harris and Zaretsky produced statistically similar 
results. 

2. For a ninth power Hertz stress-life life exponent for ball 
bearings, the Lundberg-Palmgren formula best predicts life. How- 
ever, for a 1 2th power relation reflected by modem bearing steels, 
the Zaretsky formula based on the Weibull equation is superior. 

3 . For cylinderical roller bearings, the flat roller profile with- 
out considering edge or stress concentrations has the longest 
predicted life followed by the end-tapered profile, the aerospace 
profile and the fully crowned profile, respectively. The fully 
crowned profile produces the lowest lives. While there are life 


TABLE 4 — ROLLING-ELEMENT BEARING LIFE 
ADJUSTMENT FACTORS FOR LIFE PREDICTION 
METHODS BASED ON L^Akkpwu = 1 


Confidence 

level, 

percent 

Life adjustment factor I 

Angular-contact 
ball bearings 

Deep-groove 
ball bearings 

Cylindrical roller 
bearings 

Lundberg-Palmgren method I 

90 

0.39 

0.78 

2.1 

95 

0.14 

0.47 

0.99 

99 

0.012 

0.15 

0.19 

99.9 

— 

0.03 

0.02 

Modified Ioannides-Harris method i 

90 

0.11 

0.70 

0.47 

95 

0.04 

0.44 

0.22 

99 

0.003 

0.15 

0.042 

99.9 

— 

0.04 

0.004 

! Zaretsky method 1 

90 

0.13 

0.55 

0.88 

95 

0.05 

0.33 

0.38 

99 

0.006 

0.10 

0.08 

99.9 

— 

0.02 

0.01 
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differences between the end tapered profile and the aerospace 
profile, these differences may not be significant. 

4. The effect of edge loading on the flat roller profile is to 
reduce life at high loads by as much as 98 and 82 percent at 
the low loads. The actual percentage calculated depends on the 
analysis used. 

5. The resultant predicted life at each stress condition not 
only depends on the life equation used but also on the Weibull 
slope assumed. The smallest change in predicted life with Weibull 
slope comes with the Zaretsky equation. For Weibull slopes of 

1.5 and 2, both the Lundberg-Palmgren and Ioannides -Harris 
(where r u equal 0) equations predict lower lives than the ANSI/ 
ABMA/ISOstandard. 

6. Based upon the Hertz stresses for line contact, the load- 
life exponent p of 10/3 for roller bearings results in a maximum 
Hertz stress-life exponent n equal to 6.6. This value is inconsis- 
tent with that experienced in the field. A value of n equal 8.1 for 
roller bearings should be used with the Lundberg-Palmgren 
method and 9.9 for the Zaretsky method whereby p equals 

4.05 and 4.95, respectively. 
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